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We study the nonequilibrium switching phenomenon associated with the metal-insulator transition 
under electric field E in correlated insulator by a gauge-covariant Keldysh formalism. Due to the 
feedback effect of the resistive current I, this occurs as a first-order transition with a hysteresis of 
I-V characteristics having a lower threshold electric field ( ~ 10 4 Vcm _1 ) much weaker than that 
for the Zener breakdown. It is also found that the localized mid-gap states introduced by impurities 
and defects act as hot spots across which the resonant tunneling occurs selectively, which leads to 
the conductive filamentary paths and reduces the energy cost of the switching function. 
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In correlated electronic systems, the Coulomb interac- 
tion and the electron-phonon coupling give rise to various 
long-range orderings of spin, charge, and orbital degrees 
of freedom of electrons, providing rich phase diagrams 
and intriguing phenomena such as colossal magneto- 
resistance [l[ . The collective response can be significantly 
sensitive and amplified in comparison with that in semi- 
conductors [2|, because many electrons cooperate in a 
short length scale of nanometres owing to the high elec- 
tron density. These orderings often lead to an insulating 
behavior with an energy gap in the single-electron spec- 
trum represented by the Mott gap [3|. In sharp contrast 
to the band insulators, the gap itself can be controlled 
by the external stimuli, as observed in experiments 
on the metal-insulator transition driven by the electric 
field 0, [j| or the light irradiation 0, 0|- Furthermore, 
metal-insulator switching phenomena have been observed 
in other correlated electronic systems such as organic 
charge-transfer compounds Q, La2-ajSr a; Ni04 H, one- 
dimensional Mott insulators S^CuOs/SrCuO^ [10|. Re- 
cently, the application of the switching phenomenon to 
electronic devices has also been seriously considered [ll[ . 
One important observation here is that the threshold 
electric fields observed in the correlated systems are typi- 
cally 10 4 V//cm @,G3 which is much less than ~ 10 6 V/cm 
expected from the simple Zener breakdown (see below). 
This suggests a positive feedback effect of the collective 
nature of the metal-insulator transition in the switching 
phenomena. Also the current is often non-uniform and 
confined in narrow paths or filaments 0, fl2| . 

Theoretically, on the other hand, the description of 
the nonequilibrium states still remain a challenge even 
though there are several related works [HI, [3 [l5[ . Espe- 
cially the non-perturbative treatment of the steady state 
under a strong electric field has been a difficulty. Re- 



cently, we have developed such formalism to deal with 
the far-from-equilibrium states [l6j]. This enables us to 
exactly incorporate the effects of the electric field into the 
Dyson equation for the nonequilibrium Green's function, 
which is written in a compact form by using the Moyal 
product in the gauge-covariant Wigner representation. 

In this Letter, we develop a theory of resistive switch- 
ing phenomenon in the spin/charge ordered correlated in- 
sulators employing the gauge-covariant Keldysh formal- 
ism combined with the mean-field approximation to the 
electron-electron interaction. Then, the hysteric resistive 
switching due to the applied electric field has been ob- 
tained theoretically for the first time as far as we know. 
The theory also accounts the experimentally observed 
low threshold field and filament formation. 

We study the one-dimensional interacting electrons de- 
scribed by 
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where g is assumed to be a constant and repulsive and the 
other notations are standard. This interaction naturally 
leads to the spin density wave ordering at the wavevector 
2kp (hp: Fermi wavenumber), which is assumed to be 
half of the reciprocal lattice vector G, i.e., half- filling, 
and introduces the gap 2 A. The sign of the gap is the 
opposite for the opposite spin, and we can just consider 
the two copies of the spinless electrons by the mean field 
Hamiltonian; 
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where, 
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= t (c P .R, c p.L) is the two component operator 
corresponding to the right-going and left-going electrons 
near ±kp with the dispersion ±vp (v. velocity). Here, 

2A := g(cl r c p ,l) is the self-consistently determined gap 
(see Eq. (Q}). While this ordered state in equilibrium 
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FIG. 1: The tilted band structure under an external electric 
field E. At each spatial position R, the momentum p is defined 
as shown in the orthogonal direction. The conduction band 
bottom and the valence band top cross the energy ui = at 
R2 and Ri, respectively. The localized impurity state at R = 
Rimp. with green color gives rise to the resonant tunneling. 



is well-known, we are interested in the nonequilibrium 
phase transition driven by the electric field. Note that 
the ordering is commensurate, and the phason degrees 
of freedom is quenched in sharp contrast to the slid- 
ing charge density wave problem [ijj • The interaction is 
treated in the mean-field approximation, which is justi- 
fied in the weak to intermediate-coupling regime, though 
a more elaborate treatment is required in the strong- 
coupling regime. 

We separate the problem into two steps, i.e., (i) to de- 
scribe the current flowing state under the electric field in 
the mean field Hamiltonian Eq. ((2]), and (ii) to solve the 

self-consistent equation 2A = g(Cp R c Pt L) for the gap. 
The first step is basically the Zener tunneling problem 
studied previously [H, dl HJ SH H, SI . As schemat- 
ically shown in Fig. 1, the electrons tunnel through the 
energy gap. The band structure is spatially tilted by the 
potential energy gain —eER where R is the real space 
position and — e is the electronic charge. One can con- 
sider the locally defined band structure as a function of 
the momentum (which is represented along the transverse 
axis) at each R, and the equi-energy line crosses the bot- 
tom (top) of the conduction (valence) band at R = R\ 
(R = i?2)- The wavefunctions of conduction and valence 
bands tunnel through the potential barrier between Ri 
and i?2 from the both sides. 

There are three length scales with this problem; (i) the 
correlation length £ = hv/2A associated with the energy 
gap 2A, which describes the characteristic extent of the 
wave packet relevant to the tunneling, (ii) the tunneling 
length £ = 2A/eE = R-2 — R\ over which an electron 
can gain the energy 2A by the electric field E, and (iii) 
the mean free path I. The Zener tunneling occurs quite 
differently depending on the relative magnitudes of these 
length scales. We will focus below the case of I ^> (; 



the mean free path is much longer than the correlation 
length, or the energy gap 2A is much larger than the 
energy broadening fmji due to the impurity scatterings. 
Then the Zener tunneling is controlled by the ratio of 
£/C- When £ 3> the Zener tunneling probability can 
be calculated in the semi-classical approximation as ~ 
cxp[— 7r£/2£] [23|]. As we increase the electric field so that 
£ < C the wave packet extends from i?i to i? 2 and the 
metallic conduction occurs, i.e., the crossover between 
the Zener tunneling and Ohmic regions. 

Although this picture is valid qualitatively, it is cru- 
cial to consider steady state with the dissipative cur- 
rent flowing to describe the nonequilibrium phase tran- 
sition. We perform the self-consistent calculations of 
Green's functions and self-energies in the Keldysh space 
in the gauge-covariant Wigner representation, which is 
now composed of the mechanical energy and momen- 
tum [ll| [24j]. It is necessary to introduce the Green's 
functions and the self-energies in the Keldysh space, 

n ( G R 2G< \ , „ / t R 2t< \ 

G := ( dA j, and E := ^ Q ±A j , respec- 

tively [1 61 ] . Kinetic equations of the functions are given 
in the form of the Dyson equations: (C — E) * G = 1 
and G* (£ - E) = 1 with t(u,p) := ui - H(p). The 
symbol * denotes the Moyal product: (/ 7k- g)(x) = 
{veEKvy I dydzf(y)g{z)e-^^ x>1 -y^ s ^ xV - z ^\ where 
x, y, z denote the two-dimensional energy-momentum co- 
ordinates (w, vp), f and g are the smooth functions of the 

energy-momentum, and := ( ^ ^ 

Now, we turn to the calculation of G R ' A and G < for 
the Landau-Zener model with the ^-functional random 
impurity potential. Let us start with the pure case of 
T, R ' A — ~^ir\ with an infinitesimal number r\. By using a 
Fourier transform: T\f\(u,i) := J ^E.f(uj,p)e- 2l P £ / heE , 
we obtain the Green's function in the pure case as 

$ S (T 



G R uf c (uJ,p) 



E- 

4^ eE 



dedT- 



e)$t(T-e) 2lpe 

_ £-1 IgeEft i 



T ± ir] 



(3) 



where $ s := ' are the solu- 

tions of the following Weber equations [ll|: 
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with a 



normalization condition: X) s =i 2 / d\$>(A — uj)$'(\ — 
u>') = 6((u> — u>')/eE). We employed the self-consistent 
Born approximation; T, R ' (w) = TiiU 2 g R ' A (u>,e = 0), 
with the density of impurities n, and the strength of 
the potential u leading to the lifetime r := (vh 2 / riiU 2 ) 
and the mean-free path £ := tv. Then, we calcu- 
late the retarded and advanced Green's functions 
through G R ' A = G R ± + G R ± * ± R ' A * G R ' A . Fi- 
nally, the lesser Green's function is obtained as 
G< S f F * G A - G R * f F + G R * [f F *C]*G A with 
the Fermi distribution function fp(u>). Here, we ne- 
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FIG. 2: (a): The right-hand side of the gap equation Eq. Q, 
(Cp h c p ,l)/2 in the unit of (F/lOftu) as a function of the gap 
A for C/C = eEhv/A&l = 0.0,0.8, 1.2, with Y being the en- 
ergy cutoff, i.e., the half bandwidth of the equilibrium states. 
Here, Ao is the gap in the equilibrium for g — 5hv. The solid 
straight line represents A/g in the same unit and the crossing 
of these two gives the solution(s) to the mean field equation 
((4]). (b): The obtained current / as a function of the electric 
field E. The dashed line is a guide to the eyes. The I-E char- 
acteristics clearly shows the hysteresis and switching behavior 
of the current. The red dotted curve represents the current 
obtained for the fixed gap Aq. 



gleet a vertex correction since it gives only a minor 
correction [25| . 

The self-consistent mean-field gap equation of A for a 
given interaction strength g is given by 



1 



( c p,r c p,l)/2 = 
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Ain 27th 



tr 



a x G<(u,p) , (4) 



the right-hand side of which is a function of the electric 
field E and the gap A itself. Throughout this paper, we 
take g = 5hv. The solution is obtained by the crossings 

of the straight line A/g and the curve for (c^ l c p ,r) /2 
in Fig. 2(a). There exist three regions of the strength 
of the electric field as (I) E < E u (II) E 1 < E < E 2 , 
and (m) E > E2, where the number of solutions to 
Eq. (0| is two, three, and one, respectively. Note that 
the stability of each solution is determined by the con- 
dition 9f(c^ fl Cp,i}/2j /dA < l/g. Thus, in region (I), 

the finite- A solution is the only stable one. In Fig. 2(a), 
we show the case of the equilibrium (E = 0) with the 
gap Ao. In the region (II), there are two stable solu- 
tions, i.e., A = and A ^ as shown for the case of 
C/£ = eEhv/AAl = 0.8 in Fig. 2(a), except the inter- 
mediate unstable one. This is the typical situation of 
the first-order phase transition. In the region (HI), the 
stability of the A ^ solution is lost, and the metallic 
state (A = 0) becomes the only stable solution, as shown 
for the case of £/£ = 1.2 in Fig. 2(a). Therefore, we 
conclude that the spin/charge ordered system shows the 
first-order-like switching phenomenon. 

We argue that these two threshold electric fields are es- 
sentially given by E\ = A /evr and E% — A^/ehv. It is 
easy to understand that E 2 is the Zener breakdown field 
since at E > E 2 , the gap does not prevent the metal- 



lic current flow and hence the insulating state is unsta- 
ble. To understand why the lower threshold field Ex ap- 
pears, it is useful to consider the instability of the metallic 
current-carrying state. The steady state with the current 
is characterized by the shift of the electron distribution 
function by the amount Sk — eEr with r = £/v being 
the mean-free time. With this shift, the energy differ- 
ence between the right and left-moving electrons at the 
shifted Fermi level is 5e — 2vSk. When this energy is 
larger than the gap 2A in the equilibrium state, the 
instability toward the SDW/CDW disappears. This con- 
sideration leads to the estimation E± — Ao/evr, which 
is smaller than the Zener breakdown field E% — A^/ehv 
by the factor H/(tA ) < 1. 

Now we study the physical properties associated with 
the switching phenomenon. Of the most important is the 
I-E characteristics. The current / flowing through the 
sample is obtained from the relation 
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Figure 2(b) shows the I-E characteristics corresponding 
to the first-order phase transition of the order param- 
eter obtained in Fig. 2(a). There occurs the jump of 
the current, the upper branch of which corresponds to 
the metallic conduction while the lower branch to the 
Zener tunneling in the insulating state. The two thresh- 
old electric fields E\ and E 2 can be separated by a factor 
as discussed above, and the change of the current is by 
the factor of ~ exp(4A 2 /ehvE). 

We also propose the measurement of the local density 
of states (LDOS) in terms of the scanning tunneling spec- 
troscopy (STS) to study the nonequilibrium state. Based 
on the formula for the tunneling current given by Meir 
and Wingreen [2(|, we have calculated the STS LDOS 
as shown in Fig. 3. There appears a peak at the middle 
of the gap whose height is proportional to the tunnel- 
ing current. Namely, the tunneling occurs through the 
in-gap density of states induced by the electric field. In 
the metallic state after the switching, of course the gap 
completely closes. 

Now the semi-quantitative estimation for the realistic 
situation is in order. Typically, 2A is of the order of 
leV, while H/t is ~ lOmeV^, which leads to the factor 
of ~ 100 reduction of the switching threshold from the 
Zener breakdown field ~ 10 6 V^/cm. Therefore, the ob- 
served values of the order of 10 4 F/cm[l, [l(| is in the 
reasonable range as expected from the present consid- 
eration. A threshold current density j\ is estimated by 
ji ~ (Ana)(e/S), where n denotes a density of the elec- 
tron. By using the typical values, we estimate the current 
density as j\ ~ 10 s A/cm 2 . This is a very large current 
density, and can be usually realized only in the pulse cur- 
rent experiment since the huge heat generation makes the 
sample burned out. 

However, the breakdown and switching occur often in 
the filamentary paths of metallic regions 1, G3. As 
mentioned above, the emerging in-gap state is associ- 
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FIG. 3: The local density of states (LDOS) Im [G A A j as a 

function of the normalized energy u>/ Ao (Aq: the energy gap 
in the equilibrium) for £/£=6.4 (dashed line), 3.2 (solid line) 
and 0.4 (dotted line), respectively. (£ = hv/2Ao; correlation 
length, £ = 2A /eE) 

ated with the tunneling. In real materials, there are 
often in-gap states due to the impurities, vacancies, etc 
even without the electric field. Suppose there is an im- 
purity level in the gap as shown by the green peak at 
R = -Rimp. in Fig. 1. Let t\ be the tunneling ampli- 
tude from the valence band at R = Ri to the impu- 
rity level, while t 2 be that from there to the conduc- 
tion band at R = R 2 . Note that t\ and t 2 are ex- 
ponentially small, i.e., t\ ~ exp(— 7r(i?j mp . — i?i)/2£), 
t 2 ~ exp(— 7r(i?2 — i?im P .)/2C) <SC 1, with the product t\t 2 
being the Zcncr tunneling amplitude without the reso- 
nant level. Considering the two barrier problem with 
the small tunneling amplitudes of t\ and t 2l the tunnel- 
ing probability through the two barriers has the peak 
height given by T max = [2tit 2 /(t 2 + t 2 )] 2 within the nar- 
row energy width 8e = W{t\ + t 2 ), which can be trans- 
lated into the width in the real-space 5R = Ss/(eE) in 
Fig.l. When the extent of the electron wave packet is 
larger than SR (which is the case in the limit of small 
tunneling amplitude), the averaged tunneling probabil- 
ity is of the order of t 2 t 2 /(t 2 + t 2 ). This is much larger 
than that without the impurity level, i.e., (iiia) corre- 
sponding to the tunneling amplitude t\t 2 for the Zener 
tunneling. This means that there appear "hot spots" at 
R\ and R 2 which are spatially separated by Q. There- 
fore the impurities act as the nucleation centers of the 



nonequilibrium first-order metal-insulator phase transi- 
tion discussed above. With the random configuration of 
the impurities, the current can find the path along where 
this nucleation centers populate densely compared with 
the other spatial region. This leads to the filamentary 
paths of the metallic regions as observed experimentally. 

Since the width of the filamentary path is given by 
the correlation length £, which characterizes also the 
spatial variation of the SDW/CDW order parameter, a 
threshold current I\ is estimated by 1\ ~ 7r£ 2 • j\ — 
(nav 2 1 A)(ireh) ~ 1/iA. From this expression, we found 
that the switching occurs with rather tiny current in cor- 
related insulators. Moreover, from a point of view of 
the Joule heating, the correlated insulator is more ad- 
vantageous over semiconductors. The Joule heat corre- 
sponding to the breakdown is written as P = j± ■ E\ = 
(Aa) 2 (1 / hvr) . For the correlated insulator, this is es- 
timated as P cor . ~ 10 12 VA/cm . This value is rather 
similar to that in the typical semiconductors. However, 
as mentioned above, the current is confined within fila- 
ments of nano-scale in the correlated systems, while it is 
rather uniformly distributed in semiconductors. There- 
fore, the total heat generation P^ or is expected to be 
much smaller in correlated insulator as P c ' or = Qir(, 2 j ~ 
q x 10 _2 VA/cm , where g denotes a density of the fil- 
aments. Once some filaments appear, the voltage drop 
across the sample disappears and no additional filaments 
are needed. 

To summarize, we have studied the switching phe- 
nomenon of the correlated insulator under an applied 
electric field. Due to the feedback effect of the current 
on the spin/charge ordering, the switching occurs with 
much weaker field/smaller current/ smaller heat gener- 
ation as compared with those expected from the simple 
Zener breakdown picture. This finding will be useful for 
the future application of this phenomenon to memory 
and switching devices. 
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